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Abstract: Elementary keV sterile Dirac neutrinos can be a natural ingredient of the composite neu- 
trino scenario. For a certain class of composite neutrino theories, these sterile neutrinos naturally have 
the appropriate mixing angles to be resonantly produced warm dark matter (WDM). Alternatively, 
we show these sterile neutrinos can be WDM produced by an entropy-diluted thermal freeze-out, with 
the necessary entropy production arising not from an out-of-equilibrium decay, but rather from the 
confinement of the composite neutrino sector, provided there is sufficient supercooling. 
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1 Introduction 

Sterile neutrinos with masses at the keV scale are a popular warm dark matter (WDM) candidate [1- 



22], that may potentially account for small-scale structure formation (see e.g [17, 23, 24]) and possibly 
explain large pulsar kick velocities [15, 25]. Sterile neutrino WDM can be produced non-thermally 
via (non)-resonant oscillations from the active neutrinos [2-4, 12, 14, 15, 26-29], by decays from the 
inflaton [30, 31], or thermally with subsequent entropy dilution (see e.g. [32, 33]). Typically, the 
parameter space spanned by the mass (hereafter m^) and active-sterile mixing angle (hereafter 9 c i) for 
sterile neutrino WDM is most tightly constrained by Lyman-a [27, 32] and x-ray flux [18, 26, 29, 34, 35] 
bounds, along with free-streaming, Tremaine-Gunn and big-bang nucleosynthesis bounds, too (see e.g. 
[15, 28]). The aggregate effect of these bounds depends on the production mechanism of the sterile 
neutrino WDM. In particular, at present purely non-resonant production is disfavored, while windows 
exist for resonant production, production from inflaton decay, or from entropy-diluted thermal freeze 
out [14, 15, 26, 27]. 

In this Note, we show that elementary keV Dirac sterile neutrinos can be a natural feature of the 
composite neutrino scenario [36-40], in the same way that the light fermions of the standard model 
(SM) can arise naturally in the extended technicolor framework [41]. Briefly, the composite neutrino 
scenario is a class of theories in which the right-handed neutrinos are composite bound states of a 
confining hidden sector (CHS). 

The possibility of such keV sterile neutrinos was first mentioned briefly in Ref. [42], and some of 
its x-ray flux bounds were investigated in [43]. In this Note, we present a more generalized discussion 
of this mechanism that is independent of the precise details of the confining sector, and then proceed 
to investigate the possible cosmological histories for this WDM candidate. We show certain classes 
of CHS's can naturally produce keV sterile neutrinos with active-sterile mixing angle in the reso- 
nant production window, and a freeze out temperature > TeV. Provided the post-inflation reheating 
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temperature is below the TcV scale, then these kcV sterile neutrinos could be WDM produced non- 
thermally via the usual resonant production mechanism [2-4, 12, 14, 15, 26-29], or by a combination 
of infiaton decay and subsequent non-resonant production [30, 31]. 

As mentioned above, an alternative to non-thermal WDM production is ultra-relativistic thermal 
production followed by entropy dilution (see e.g. [32]). This has the advantage of producing colder 
WDM than resonant production and can better evade the Lyman-a bounds. Usually the diluting 
entropy is produced by the out-of-equilibrium decay of a sufficiently long-lived heavy particle. In 
this Note we examine another compelling possibility: The first-order phase transition induced by the 
confinement of the hidden sector can also produce significant entropy if there is sufficient supercooling. 
This results in thermal keV WDM. We will discuss the details of this mechanism. 

2 The Composite Dirac Neutrino Model 
2.1 Setup 

The generic theory of interest is a low-energy effective field theory below a scale M. Its group structure 
is G c ® Gf ® Gsm, with G c a confining group called z^-color, Gsm the SM gauge groups (or a UV 
extension), and Gp a global (or weakly gauged) hidden flavor group. The theory consists of three 
sectors 

X~G C ®G F , £~G P , q~G SM ®G F , (2.1) 

and which interact only via M-scale irrelevant operators. We call x 'preons' and say they belong to 
the CHS. Here q denote the SM fields extended to also carry hidden flavor Gf, and we say £ comprise 
the 'extended hidden sector' (EHS). We assume that the \ and £ are purely chiral fermions, but we 
emphasise that like the SM sector, the x ancl £ ma y consist of various different irreps. 

The i/-color group confines at a confinement scale A -C M. Necessarily M 3> v, the electroweak 
scale, so it is convenient to define two parameters 

e = A/M<l, 6 = v/M<l. (2.2) 

Confinement of the CHS produces preonic bound states, which we shall crudely denote as x p '- The 
superscript denotes the number of preons participating in the bound state. Formation of a scalar 
condensate x m with (x m ) 7^ generically induces a spontaneous breaking of the hidden flavor group 
Gf — > G' F C Gf- This produces a new sub- A effective field theory, which consists of: preonic bound 
states; £ and q decomposed into G F irreps; and also light 'hidden pions'. There are three crucial ideas: 

(i) If the CHS has non-trivial G' F anomalies, then anomaly matching of the CHS to its confined 
phase, with £ and q acting as chiral spectators, implies that there are massless fermionic bound 
states after confinement. The remaining bound states generically have masses ~ A, except for the 
hidden pions, which can be massless or have arbitrarily small masses, depending on the nature of the 
Gf symmetry breaking. We assume the pion masses are sufficiently small that they make negligible 
contributions to the DM energy fraction. 

Hereafter we shall assume G' F = U(1)f, and that there are precisely three massless bound states 
all with the same U(1)f charge 1 . For simplicity we assume the massless bound states have the same 
number of preons, hereafter denoted n, necessarily an odd integer. We shall suggestively denote these 
bound states as njj, i = 1, 2, 3 with U(l)p charge F{np) — +1. Explicit examples of preonic theories 

1 In this case decomposition of q under Gp —> U(l)p could result in multiple copies of SM irreps, also with the same 
U(1)f charges, which could be the source of flavor. 
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capable of producing such spectra are presented in Ref . [42] . The corresponding sub- A EFT that we 
shall consider hereafter is shown in Table 1. In producing this EFT, we require that the mechanisms 
of Gf — > U(1)f breaking and electroweak symmetry breaking are independent, at least to a good 
approximation. 
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Table 1. U(1)f charges assignments to the massless bound states n_R and the SM fields q = {<j>, Q, U, D, L, E}, 
which also have the usual SM charges (not shown). The tir are SM sterile by construction. 

One can check 2Y — F = B — L, so U(l)p is nonanomalous, and the electroweak symmetry breaking 
(EWSB) pattern is 

SU(2) L ® U(l)y ® U(1) F -> U(1) EM ® V(1) B _ L . (2.3) 

That is, one obtains Dirac neutrinos, with the tir acting as right-handed neutrinos. Note U(1)f may 
be gauged, but we assume its gauge coupling and kinetic mixing with the photon are sufficiently small 
that they can be neglected. 

(ii) For the sub- A EFT in Table 1, there exist irrelevant operators that couple the preons of the 
massless ur - i.e. the G c singlets x™ ~ to the SM singlet L^tp. Such an operator is generically of form 

M^m l ^X n ^-^L^ur , (2.4) 

after confinement. That is, this operator produces a suppressed Yukawa in the sub-A EFT. Since tir 
are massless and there is B — L symmetry (2.3), this operator leads to light Dirac neutrino masses 
after EWSB, compared to the electroweak scale. 

There may also be other vector-like right-handed fermionic bound states Nr and iV£, with 
F(Nr,l) = +1 We shall again assume for simplicity they contain n preons. Such bound states 
must form Dirac fermions with A scale masses, and the Nr will generically also have operators of form 
(2.4). Nr^l are therefore A-scale sterile Dirac neutrinos. 

(iii) Under decomposition into U(l)p irreps, the chiral EHS fields £ may form real U(l)p represen- 
tations and acquire masses. However, because the EHS couples only irrelevantly to the condensate vev 
(x m ) responsible for Gf — > U(1)f, the mass terms must be suppressed. This is the same mechanism 
which suppresses the quark and lepton masses in Extended Technicolor theories [41]. Explicitly, for a 
Dirac fermion £,r^l, such mass terms arise from operators of the form 

_L_^ Ke^'HUn , (2.5) 

after confinement 2 . If also F(^r^) — +1, then there may exist irrelevant operators that couple the 
corresponding G c singlet x m £ to L^cp, noting any renormalizable coupling of £ directly to L^cj) is 
forbidden by the Gf chiral structure. That is, we could have 

^ W>£ -> e^L L kn ■ (2.6) 



2 There may also be mass cross terms involving for example. However, we assume that such cross-terms, i.e 

involving composite and elementary states, are suppressed by the details of the UV theory above M. An analogous 
assumption must also be made for the proton decay operator uude/M 2 . 
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Consequently, such a £r,l forms an elementary sterile Dirac neutrino with naturally suppressed mass 
term ~ Ae^ 3m ~ 2 ^ 2 and coupling to the active sector ~ e 3m / 2 . j n p r i nc ipl e , there may be several 
species of such a Dirac neutrino, as well as other EHS fermions with F ^ ±1 that acquire Dirac or 
even Majorana masses of the same size. 

2.2 Spectrum 

We may classify the sub- A EFT by a tuple (n, to), where n (odd > 3) is the number of preons in the 
sterile neutrino bound states, and to (even > 2) is the number of preons in the symmetry breaking 
condensate. After EWSB, from eqs. (2.4)-(2.6) a (n,m) theory has neutrino mass term, 



(2.7) 



where vl is the SM active neutrino. Each entry of this mass matrix denotes the prefactor of an 0(1) 
sub-block, whose dimensions depends on the number of species of each type of sterile neutrino. For 
example, the upper left entry must be 3 x 3. 

For to < n — 1, the mass spectrum can be determined by expansions in e and 9. One obtains at 
leading order 

to/ ~ ve 5 , m d ~ Ae^~ , rrih ~ A . (2.8) 

Here the superscripts I, d and h denote 'light', 'dark' and 'heavy'. The left-handed mass basis is, at 
leading order in e and 9, 








n 
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(2.9) 



and the right-handed mass basis is 



/ I 

2, 



3(5 

"e 



\ e 2 e 3n ~ 5 e 2 e 3 -^^- 1 

We emphasise that eqs. (2.9) and (2.10) denote only sub-block prefactors; the entries of the sub-blocks 
themselves are generically 0(1) numbers multiplied by the appropriate prefactor. 

It is clear from eq. (2.9) that the dark-active mixing angle 9d ~ 9. One can then rearrange eqs. 
(2.8) and (2.9) into 




1 02 £ ~^ Ufl ■ (2-10) 



TO/ \ A /to/ 



m d e d ~ v -!• , — , (2.11) 

V v ) m d \ v J 

in which the right-hand sides are fully specified by (n, m) and the requirement that to/ ~ 0.05 eV, 
v ~ 174 GeV. Figure 1 shows sin 2 (29d) up to 0(1) uncertainty as a function of m d , with m = n — 1. 
Theories with to < n — 1 have much larger mixing angles, and are therefore ruled out by x-ray flux 
constraints, so we consider only (n, n— 1) theories henceforth. For such theories M ~ 2x 10 4 (to £ //5 keV) 
TeV, and we provide the corresponding A and e in Table 2. 

It is amusing to note that for the (n,n — 1) theories m d ~ 5 keV implies sm 2 {29 d ) ~ 3 x 10~ 10 , 
which matches the (as yet unconfirmed) Chandra results in the Willman I dwarf galaxy [45] . 
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Figure 1. Mixing angle sin 2 (2^) up to 0(1) uncertainty (light gray) as a function of mj, for (n,n — 1) 
theories. Also shown: Non-resonant production contours (dashed lines), labelled by the ratio of v d and DM 
energy fractions, ild/fiDM [9, 14, 15]; resonant total DM production contours (dash-dotted lines) for lepton 
asymmetries Yal = 8, 12, 16, 25 x 10 -6 (resp. top to bottom), and their corresponding Lyman-a lower bounds 
on the WDM mass (black dots) [14]; the Lyman-a exclusion for thermally produced WDM with subsequent 
entropy dilution (hatched region, see e.g. [15, 27] and eq. (3.7) below) assuming 100% v d WDM; the x-ray 
flux exclusion for 100% v d WDM fitted from most stringent archival data (heavy black line, see e.g. [14, 26]) 
and from the most recent observations of dwarf spheriodal galaxies [44] (heavy broken line). 
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Table 2. Confinement scale A and e for (n, n — 1) theories. Such theories with n > 7 have e 1, and are not 
considered further. 

2.3 Dirac vs Major ana 

The keV sterile neutrinos in this Note are Dirac, in contrast with the Majorana sterile neutrinos 
often considered in other WDM scenarios. The WDM production mechanisms that we consider below 
produce dominantly symmetric DM - the resonant production mechanism requires an asymmetry in 
the proper number density (n v — n j? )/n v < 1CP 2 [4, 12] - so that the DM particles and antiparticles are 
present in the same abundances to a very good approximation. The x-ray flux bounds due to sterile 
neutrinos are therefore insensitive to the mass structure, since decay modes to the active neutrino and 
antineutrino are present in both cases: I.e, the x-ray flux is due to either N — > 1/7 and N — > v c "f for 
a Majorana neutrino N, or v d — > i/j and D d — > for the present scenario. Similarly, (non)-resonant 
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production by conversion from the left-handed active neutrinos will produce the same sterile neutrino 
energy fraction, fi^, regardless of the Dirac or Majorana nature of the masses. In Fig. 1 we therefore 
use the existing results for both the x-ray bounds and production processes, without any alteration 
for the Dirac mass structure. 

The x-ray bounds could also be altered by exotic v d — > Xj decay channels, that might arise 
from M-scale irrelevant operators. We emphasize that the chiral and composite structure of the 
composite neutrino framework ensures any such operators are of sufficiently high dimension that 
the corresponding decay rates are negligible. For example, in the (3, 2) theory v d — > jv l or v d — > 
could also arise from ^X^F^a^^/M 7 which confines to (A 5 /M 7 )n fl (A + U)F^a^^ L . This 
respectively produces decay rates ~ e 12 m d /M 2 or ~ e 10 m|/M , that are negligible compared to the 
decay through mixing with the active neutrinos. 

2.4 Decoupling 

Our knowledge of the generic structure of the non-renormalizable operators permits us to consider the 
cosmological histories of the CHS and EHS, and therefore determine whether the v d sterile neutrinos 
can be a WDM candidate: satisfying the (mj, Od) bounds is necessary but not sufficient for this. For the 
(n,n — 1) theories, we now enumerate various important processes and their freeze out temperatures, 
Tf r . We assume the effective degrees of freedom at the TeV scale ~ 10 2 . 

(i) XX O YY, where X,Y € {q^,x}- These processes couple the SM, CHS and EHS. The 
dimension-5 operator <jy<j)XX is heavily suppressed, since X are all chiral. The leading operators are 
then the dimcnsion-6 



1 

M 2 



X^XY lp Y ; T fr 



^ - TeV , (2.12) 



and similarly for (fyd^cpX^X/M 2 . Note that the current collider constraint on the dark matter - 
quark interaction is insensitive to the coupling due to the large mediator mass, M [46, 47]. 

(ii) 211, where n denotes the hidden pions. This process is generated by the non-linear 

sigma operator 

gl A 4 



m d ^ L e m ' A ; T fr 



TeV , (2.13) 



{m d ) 2 M pl _ 

for the (3,2) theory, and much larger for (5,4) and (7,6). 

(iii) VjVj^ 4-> qq. This can occur also through W and Z exchange, and must freeze out before 
BBN. The pertinent operators are 



g(0d) 2 d + d gOd dTirp . rp 
2c w ' ' 1 



(0 d ) 4 M pl 



-1 1/3 

TeV 



(iv) v l L v l R O 2n. This must also freeze out before the BBN epoch. The non-linear sigma coupling 
of v L R to the hidden pions is suppressed by both the left and right mixing between active and sterile 
sectors. From eqs. (2.9) and (2.10) this leads to an extra prefactor of {Od} 3 for the non-linear sigma 
operator in eq. (2.13), and therefore a decoupling much larger than the TeV scale. 

(v) 2h «-> 2n. This is generated by the operator $ 4>{x m )^ X m I M Zm ~ 2 which confines to the 
dimension-4 operator e 3m_2 ^t^nn. This becomes efficient only below 



T fr - e 6m - 4 M p i/ 5 | < 1(T 7 eV , (2.14) 
for (n,n — 1) theories, and therefore does not produce significant recoupling. 
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3 Warm Dark Matter 



3.1 Non-Thermal WDM 

The moral of the above analysis is that approximately below the TeV scale, the SM, CHS and EHS are 
decoupled. From Table. 2, confinement of the CHS also occurs at latest at the TeV scale. As a result, 
we may imagine a scenario in which the post-inflation reheating temperature T r h < TeV. In this case, 
the sterile Dirac neutrinos v d might never be in thermal contact with the SM plasma, and therefore 
be produced non-thermally through the (non)-resonant production mechanism [2-4, 12], forming the 
WDM. 

As can be seen in Fig. 1, the predicted (rrid, 9d) values fall outside the fid > ^dm/2 non-resonant 
production region, which itself is ruled out by the combination of Lyman-a [27, 48] and x-ray flux 
bounds [14, 26]. However the {rrid,0d) ranges still overlap an allowed window for full WDM resonant 
production if there is a sufficiently large lepton asymmetry [4, 11, 12, 14, 15]. Alternatively, in this 
low reheat scenario, coupling of the sterile neutrinos to the inflaton - an SM singlet - could result in 
significant non-thermal WDM production from its decay [30, 31], with the remaining fraction (if any) 
produced by non-resonant production. 

Just as for exotic x-ray decay channels, the chiral structure of the SM, CHS and EHS generically 
suppresses the operators that may produce non-thermal sterile neutrino WDM from SM decays. For 
example, from eqs. (2.6) and (2.8) it is clear that Higgs to sterile neutrino decay rate is suppressed 
by e 3m ~ (mi/v) 2 for (n,n — 1) theories, so there is no significant production from the Higgs decay 
channel. Along similar lines to the inflaton scenario, one might putatively extend the Higgs sector 
with a SM singlet that can decay to the EHS without such suppression (see e.g. [49, 50]), however we 
do not make any such assumptions about the SM Higgs sector here. 

One might also consider production via lepton or hadron decays such as r — > e££ or B — > K££ 
respectively. The chiral structure ensures such processes can only be mediated by operators of the 
form 

^-VYfU ■ (3.i) 

This type of operator necessarily produces FCNCs, too, but the large mediator scale M easily evades 
the present bounds for quark FCNCs [51]. One finds for the dominant top decay process T/H(rrit) < 
10~ 4 . For semi-relativistic tops in thermal equilibrium, this produces a sterile neutrino energy fraction 
fid ~ 1% Qdm, so that this production channel can be neglected. Similarly production from spin-1 
bound state decays like po — > is negligible due to suppression of the rate by a (A qcc i/M) 4 factor. 

3.2 Thermal WDM 

The (3, 2) theory exhibits the interesting feature that the decoupling temperature of the EHS, T^, 
the confinement temperature of the CHS, T c ~ A, and decoupling of temperature the CHS, T x , all 
occur at the TeV scale. In contrast to the non-thermal resonant scenario, for a (3, 2) theory one may 
plausibly consider a scenario in which all three sectors are initially in thermodynamic equilibrium, the 
lepton asymmetry is small, and 

T d >T c >T x . (3.2) 

In this scenario, the EHS fermions £ freeze-out ultra-relativisitically before confinement, and there 
is no subsequent resonant production: from Fig. 1 we see that fractional non-resonant production 
at the 10% f^DM level may still occur, but we shall neglect this henceforth as it is a subdominant 
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contribution. Defining Y = n/s - the ratio of the comoving number density and entropy density - 
then for each Dirac £ species 

y, = • <«) 

where gf s is entropic effective equilibrium number of degrees of freedom at freeze-out. 

Even if only one species of £ - the Dirac £_r,l ~~ obtains a mass mj, which we assume henceforth, 
such a Y% leads to over-closure unless g^ s ~ 10 4 . This is unnaturally large since g* s ~ 10 2 for the 
SM at this scale. However, if after freeze-out the entropy increases by a factor 7, then the frozen 
out species are diluted, Y% — > The present-day energy fraction for the Dirac i/ d , which are an 

admixture dominantly composed of is then 



Y ( m d s a _ 1.1 x 10 4 / m d 



(3.4) 



in which we used So — 2.89 x 10 3 cm 3 , p c ~ 10.5ft 2 cm 3 keV, and Odm = 0.105ft- 2 . It is clear that 
we need <7* s 7 > 10 4 for a DM candidate. 

3.3 Supercooled Confinement 

The ordering (3.2) permits us to consider the confinement of the CHS as the source of entropy that 
dilutes Y% after freeze-out. The entropy production from a confinement-induced first-order phase 
transition can be significant if it occurs suddenly after supercooling [52, 53]. That is, if the confinement 
phase transition (CPT) begins at a cooler temperature Tj < T c , and the duration of the transition 
r c -c l/.ff(Tj), the Hubble time at temperature T t . 

Before confinement - at temperature T, - and after confinement - at temperature Tt > T x -, we 
suppose that we have equilibrium plasmas. By construction 

g* s (T i ) = gi s =g™+g° s ^2xW 2 , 

g« s (T f ) = g£ s =g™+g%~10 2 . (3.5) 

Here gf^ 1 , gl s and g^ denote the effective equilibrium relativistic degrees of freedom in the SM, CHS 
and the bound states. By construction, for three ur we have g^ s — 2 • 3 • (7/8) + iVn with iVn the 
number of hidden pions. We have assumed g£| ~ 10 and gf™, g% s — 10 2 . Note that since the frozen 
out £l,r have only four degrees of freedom, then g^ s ~ gl s . 

Since Tf > T x , then such entropy production leads to reheating of both the CHS and SM, because 
they only decouple later at T x . This mutual reheating means the present DM temperature, Tj, 
compared to that of the active neutrinos, TjJ?, is just 

T° (a* n v \ 1/3 / in 75 \ 1/3 

/ g* s y* s \ _ / lu 1 /g 



15 V 7.< s <4 M / V 1 ■ 1 x 10 4 (m d /5 keV) 

from eq. (3.4) and since g{ s ~ 5™. Equation (3.6) implies the entropy-diluted thermal WDM 
is red-shifted compared to the active neutrino plasma. The Lyman-a bounds [27, 32, 48] require 
non-resonantly produced WDM - at present temperature T® - to satisfy m nrp > 10 keV. Since the 
free-streaming length Aps oc T/m (see e.g. [15]), this Lyman-a bound translates to m d > 10(Xj/T°) 
keV. Together with eq. (3.6) we find that thermally produced v d may safely avoid the Lyman-a bound, 
provided 

m d > 1.5 keV . (3.7) 
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This is the Lyman-a bound displayed in Fig. 1. 

Note also that the hr and hidden pion contribution to the effective number of neutrino degrees of 
freedom, 5N° S , at the big-bang nucleosynthesis (BBN) epoch is 

SNf = (8/14)< s (g»Jg s ™) 4/3 < 0.26« s /10) • (3.8) 

It is amusing to note that the right-handed neutrinos together with the hidden pions can supply 
sufficient effective degrees of freedom at the BBN epoch to significantly contribute to the observed 
5N° S ~ 1 excess (see e.g [54, 55]). In contrast, this is difficult to achieve with seesaw models, or even 
ad hoc Dirac neutrino models. 

3.4 Entropy Production Estimate 

The massive bound states typically have masses xA, with x > 1, so they are non-relativistic. Their 
corresponding widths are generically also r ~ A. This leads to T/H{Ti) ~ M p \A/Tf ^> 1. In 
contrast, the longest-lived heavy bound state we could contemplate decays only via exchange of an 
M-scale boson, like the electroweak decay of the A baryon of QCD. In this case, the decay rate is 
r - Ax 5 e 4 . For the (3,2) theory e - 10~ 4 , so that T/H(Ti) > x 5 e 4 M pl A/T? > 1. This means that 
even for a sudden CPT, the heavy bounds states all decay within t c and generically, predominantly 
produce hidden pions and ur with energies ~ T c . It seems reasonable, then, to treat the CPT 
as a quasiequilibrium process, in which the non-relativistic heavy bound states have exponentially 
suppressed number and energy densities, while pions and tir are thermal with temperature T c . 

With this in mind, one can estimate the amount of entropy production by treating the CPT as 
a first-order phase transition in (7* s , as a function of ( = (RT) 3 . Here R is the universe scale factor 
and T the equilibrium temperature. The picture is that confinement begins at supercooled plasma 
temperature Tj, and suddenly produces the relativistic pions and riR at temperature T c , so that g* s 
undergoes a jump at Q — (R{Ti) 3 from gl s to 

g j .' //: x: • g::(r, '/,)''. (3.9) 

This expression for g£' s follows just from the definition g* s (T) = ^ Q g" (T a /T) 3 , a sum over species at 
different temperatures. After the phase transition, the plasma undergoes an adiabatic thermalization 
?*s = 9*s an d T — Tf. SM-CHS decoupling at T x 

SM / bs 

* s /y* s 



until g its = gl s and T — Tf. SM-CHS decoupling at T x follows thereafter. Figure 2 shows this history. 
Provided (T c /Ti) 3 ^> g^f/g^ ~ 10, the entropy production estimate from eq. (3.9) is then 



th = &„&(Tc\ 3 . (3 . 10) 

The important feature of this naive estimate is the (T c /Ti) 3 dependence of the entropy production. A 
more careful treatment in Ref. [52] produces the result 

One also finds Tf — [(r — 1)/3] 1 / 4 T C . Using this result and eq. (3.4), and fixing r = 2, it follows that 
for firf < firjM (i.e. "fgt s > 1.1 x 10 4 m d /5 keV) we require 

T c [2 x 10 2 \ 1/3 / m d \ 1/3 , , 

5T i0 (— ) (mSv) ■ (3 - 12) 



- 9 - 




Figure 2. A sketch of the thermal history. Species freeze-out (a-b) along the Si adiabat (lower dashed), 
is followed by the CPT (6-c), which is a first-order <?„ s phase transition in £. The CPT is followed by 
thermalization (c-d) along the S; adiabat (upper dashed) until g, s = g( s at which T — Tf. Once T = T x , the 
CHS and SM decouple. 



Note Tf — 0.76T C here, so it is plausible that Tt > T x . By comparison to eq. (3.12), the QCD maximal 
supercooling is T c /Tj ~ 1.7 [52]. However, given that this upper bound will be sensistive e.g. to the 
tunneling probabilities between the metastable (Gf symmetric) and stable (G F symmetric) vacua, the 
degree of supercooling required in this estimate is not implausible. 

4 Conclusions 

Within the composite neutrino framework, we have shown in this Note that keV sterile Dirac neutri- 
nos can be naturally produced with mixing angles appropriate for non-thermal resonant production, 
provided the composite neutrinos are all comprised of n preons and the scalar condensate vev has 
n — 1 of them. Alternatively, for a (3, 2) theory, a single keV sterile Dirac neutrino species could be 
WDM produced by entropy-diluted ultrarelativistic freeze-out. In this latter case the entropy can be 
provided by a supercooled confinement-induced phase transition. 
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